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ABSTRACT 

yThe  motion  of  a  two-dimensional  bubble  rising  at  a  constant  velocity  U 
in  a  tube  of  width  h  is  considered.  The  acceleration  of  gravity  g  and  the 
surface  tension  T  are  included  in  the  dynamic  boundary  condition.  For  T  * 
0,  the  results  of  Garabedian  and  Vanden-Broeck  show  that  a  solution  exists 
for  all  values  of  the  Froude  number  F  *  U/(ghl^^  smaller  than  a  critical 

-r  ,■  ,  .  v.  wtl  jt~L  to 

value  F  \  0.36.  In  this  paper  accurate  numerical  solutions  with  T  ¥  0  are 
c 

computed  by  series  truncations.  The  results  show  that  for  each  value  of 

. r  f  i  •  ' 

T  f  0,  there  exists  a  countably  infinite  number  of  solutions.  Each  of  these 
solutions  corresponds  to  a  different  value  of  F.  As  T  tends  to  zero,  all 
these  solutions  approach  a  unique  limiting  solution  characterized  by 
F  =  F*  0.23.  Therefore,  the  degenerancy  associated  with  T  *  0  is  removed 
by  including  the  effect  of  surface  tension.  In  addition  the  profile  corre¬ 
sponding  to  F  =  F*  is  found  to  be  in  good  agreement  with  experimental  data. 
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SIGNIFICANCE  AND  EXPLANATION 


The  motion  of  a  two-dimensional  bubble  rising  at  a  constant  velocity  U 

in  a  tube  of  width  h  is  considered  (see  Figure  1).  The  acceleration  of  the 

gravity  g  and  the  surface  tension  T  are  included  in  the  formulation  of  the 

problem.  For  T  ■  0,  Garabedian  and  Vanden-Broeck  have  shown  that  a 

solution  exists  for  all  values  of  F  smaller  than  a  critical  value  Fc. 

Garabedian1  used  a  stability  argument  to  suggest  that  the  only  physically 

2 

acceptable  solution  is  the  one  corresponding  to  F  *  Fc*  Vanden-Broeck 
solved  the  problem  with  T  *  0  numerically.  He  found  that  Fc  *  0.36.  This 
value  is  about  forty  percent  higher  than  the  experimental  value  of  0.25  of 
Collins3* 

In  the  present  paper  we  show  that  the  discrepancy  between  the  theoretical 
value  0.36  and  the  experimental  value  0.25  is  considerably  reduced  by  taking 
into  account  the  effect  of  surface  tension .  We  show  that  for  each  value  of 
T  *  0,  there  exists  a  countably  infinite  number  of  solutions.  As  T  tends 
to  zero,  all  these  solutions  approach  a  unique  limiting  configuration 
characterized  by  F  «  F*  ~  0.23. 

These  results  imply  that  the  physically  revelant  solution,  when  T  *  0, 
is  not  the  solution  corresponding  to  F  *  Fc  -  0.36  but  the  soltuion 
corresponding  to  f  *  F*  -  0.23. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  author  of  this  report. 


RISING  BUBBLES  IN  A  TWO-DIMENSIONAL  TUBE  WITH  SURFACE  TENSION 


Jean-Marc  Vanden-Broeck 


I.  Introduction 

We  consider  a  bubble  rising  at  a  constant  velocity  U  in  a  two- 
dimensional  tube  of  width  h.  We  choose  a  frame  of  reference  moving  with  the 
bubble  and  we  assume  that  the  bubble  extends  downwards  without  limit  ( see  Fig. 
1.) 

As  we  shall  see  the  problem  is  characterized  by  the  Froude  number 


and  the  Weber  number 


U 


(gh> 


1/2 


pU2h 


(1) 

(2) 


Here  g  is  the  acceleration  of  gravity,  T  the  surface  tension  and  p  the 
density  of  the  fluid. 

Garabedian  [1]  considered  the  problem  with  T  ■  0  (i.e.  a  ■  «).  He 
presented  analytical  evidence  that  a  solution  exists  for  all  values  of  F 
smaller  than  a  critical  value  Fc.  He  then  used  an  energy  argument  to  suggest 
that  the  only  physically  significant  solution  is  the  one  for  which  F  =  Fc. 

In  addition  he  showed  that  Fc  >  0.2363  and  guessed  the  value  Fc  3  0.24. 
Birkhoff  and  Carter  [2]  and  Collins  [3]  using  different  approaches  obtained 
the  value  0.23. 

Vanden-Broeck  [4]  solved  the  problem  with  T  «  0  numerically.  His 
results  confirm  that  a  solution  exists  for  all  values  of  F  smaller  than  a 
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Sketch  of  the  flow  end  the  coordinates.  The  broken  lines  correspond  to  the 
reflexion  of  the  flaw  domain  into  the  wall  y  *  ^  . 


critical  valua  ?c .  However  Vanden-Broeck ' a  results  do  not  support  the 
approximate  calculations  of  Garabedian  [1],  Birkhoff  and  Carter  [2]  and 
Collines  [3].  Vanden-Broeck  found  Fc  -  o.36.  This  value  is  about  forty 
percent  higher  than  the  experimental  value  0.25  presented  by  Collins  [3] . 

In  the  present  paper  we  show  that  the  discrepancy  between  the 

experiemntal  value  F  -  0.25  and  the  theoretical  value  Fc  “  0.36  is 
considerably  reduced  by  taking  into  account  the  effect  of  surface  tension.  We 
show  that  for  each  value  of  T  *  0,  there  exists  a  countably  infinite  number 
of  solutions.  Each  of  these  solutions  corresponds  to  a  different  value  of 
F.  As  T  tends  to  zero,  all  these  solutions  approach  a  unique  limiting 
configuration  characterized  by  F  ■  F*  ~  0.23.  The  corresponding  profile  is 
found  to  be  in  good  agreement  with  the  experimental  data  of  Collins  [3] . 

It  is  interesting  to  note  that  the  present  inviscid  problem  is 
qualitatively  similar  to  the  problem  of  viscous  flow  in  a  Hele-Shaw  cell 
(Saffman  and  Taylor  C5],  McLean  and  Saffman  [6],  Vanden-Broeck  [7]).  Both 
problems  are  characterized  by  a  continuum  of  solutions  for  T  -  0  and  a 
discrete  set  of  solutions  for  T  >  0. 

The  problem  is  formulated  in  Sec.  II.  The  numerical  procedure  is 
described  in  Sec.  Ill  and  the  results  are  discussed  in  Sec.  IV. 


II.  Formulation 


Let  us  consider  the  steady  two-dimensional  potential  flow  of  an  inviscid 
incompressible  fluid  past  a  bubble  in  a  tube  of  width  h  ( see  Fig.  1 ) .  The 
pressure  in  the  bubble  is  assumed  to  be  constant.  We  introduce  cartesian 
coordinates  with  the  origin  at  the  top  of  the  bubble  and  we  assume  that  the 
bubble  is  symmetric  about  the  x-axis.  Gravity  acts  in  the  negative  x- 
direction.  As  x  ♦  •,  the  velocity  approaches  the  constant  U. 

We  define  dimensionless  variables  by  taking  U  as  the  unit  velocity 
and  h  as  the  unit  length.  We  denote  the  potential  function  by  $  and  the 
stream  function  by  i|>.  In  addition  we  introduce  the  complex  velocity  by 
5  *  u  -  iv  and  we  define  the  function  x  -  10  by  the  relation 

C  «  u  -  iv  *  eT”^®  (3) 

Without  loss  of  generality  we  choose  $  •»  0  at  x  **  y  =  0  and  i|>  *  0 
on  the  surface  of  the  bubble.  It  follows  from  the  choice  of  the  dimensionless 
variables  that  t|i  =  -  on  the  wall  y  *  -j. 

we  satisfy  the  kinemetic  condition  on  the  wall  y  *  -j,  by  reflecting  the 
flow  in  that  boundary  (see  Fig.  1).  We  shall  seek  the  complex  function 
x  -  i0  as  an  analytic  function  of  f  =  $  +  i\|>  in  the  strip  -1  <  iji  <  0.  The 
complex  plane  is  sketched  in  Fig.  2. 

On  the  surface  of  the  bubble,  the  Bernoulli  equation  yields 

1/2  q2  +  gx  -  =  B  on  SJ  and  S'J’ .  (4) 

Here  q  is  the  flow  speed,  K  the  curvature  of  the  free  surface  and  B  the 
Bernoulli  constant.  In  dimensionless  variables  (4)  becomes 


Here  F  and  a  are  the  Froude  number  and  the  Weber  number  defined  by  ( 1 )  and 
(2)  respectively. 

It  is  convenient  to  eliminate  x  and  B  from  (5)  by  differentiating  (5) 
with  respect  to  $ .  Using  the  relation 


1*  +  ih.  =  _J _ -  e-T+i9 

3$  3$  u-iv 


we  obtain 


e~  t$  +  ~ T  tcoS  6  '  a  3^*T  3*)  "  0  on  **  and  S'J’ 

F 


The  kinemetic  conditions  on  IS  and  IS'  yield 


(6) 

(7) 


0-0,  4>  -  -1 

4>  <  o 

(8) 

0*o,  4*  *  0 

4>  <  o. 

(9) 

Finally  the  flow  configuration  of  Fig.  1  is  characterized  by  a  stagnation 
point  and  a  horizontal  slope  at  the  top  of  the  bubble.  This  yields  the 
additional  conditions 


'--'•-7 

at 

♦  *0,4»* 

0 

(10) 

T  -  0  -  ~ 

at 

^  *  0,  \(»  = 

-1. 

(11) 

This  completes  the  formulation  of 

the 

problem  of 

determining 

T  -  i0 

This  function  must  he  analytic  in  the  strip  -1  <  41  <  0  and  satisfy  the 


conditions  ( 7 ) - ( 1 1 ) 


III.  Numerical  Procedure 


We  solve  the  problem  by  using  the  procedure  derived  by  Vanden-Broeck  [7] 
to  Investigate  the  effect  of  surface  tension  on  the  shape  of  fingers  In  a 
Hele-Shaw  cell.  Thus  we  define  a  modified  problem  by  replacing  (10)  and  (11) 
by 


T  - 

T  - 


0  -  Y 

at  $  ■  0, 

<|>  -  0 

(12) 

0  -  -Y 

at  $  -  0, 

-  -1. 

(13) 

Here  y  Is  a  parameter  to  be  found  as  part  of  the  solution. 

We  will  solve  the  modified  problem  defined  by  (7)-(9),  (12)  and  (13),  and 
obtain  solutions  for  all  values  of  F  and  a.  Then  we  will  obtain  the 
solutions  of  the  original  problem  by  selecting  among  the  solutions  of  the 
modified  problem  those  for  which  Y  **  y. 

Following  Birkhoff  and  Carter  [2]  and  Vanden-Broeck  [7J  we  define  the  new 
variable  t  by  the  relation 


e 


-irf 


(14) 


This  transformation  maps  the  flow  domain  onto  the  unit  circle  in  the  complex 
t-plane  (see  Fig.  3). 

We  note  that 

1/3 

C  ~  (£n(1  *  it)]  as  t  ♦  *  1 

Z  ~  (1  *  t)^  as  t  +  +  1 

(see  Birkhoff  and  Carter  [2]  for  details).  Therefore  we  define  the  function 
ft(t)  by  the  relation 


-7- 


5  =  -[-An  c(l+t2)]  1/3(-£n  c)"1/3(  1-t2)Y  2/ir  [1+H(t)  ]  .  (15) 

Here  C  is  an  arbitrary  constant  between  0  and  0.5.  We  choose  C  ■  0.2.  The 
function  £2(t)  is  bounded  and  continuous  on  the  unit  circle  and  analytic  in 
the  interior.  The  conditions  (8)  and  (9)  show  that  Ci(t)  can  be  expanded  in 
the  form  of  a  Taylor  expansion  in  even  powers  of  t.  Hence 

eT~i9  =  -[An  c( 1+t2 ) 1/3 (-An  c)"1/3( 1-t2)Y  2/*(l+  \  a  t2”)  (16) 

n=1  n 

The  function  (16)  satisfy  (8),  (9),  (12)  and  (13).  The  coefficients  an  and 
the  parameter  y  have  to  be  determined  to  satisfy  (7)  on  SJ.  The  condition 
(7)  on  S'J1  will  then  be  automatically  satisfied  by  symmetry. 

We  use  the  notation  t  «  |t)ei0  so  that  points  on  SJ  are  given  by 
t  =  e10,  0  <  a  <  y.  Using  (14)  we  rewrite  (7)  in  the  form 

dS 
da 

Here  t(o)  and  6(a)  denote  the  values  of  x  and  6  on  the  free  surface 
SJ. 

We  solve  the  problem  approximately  by  truncating  the  infinite  series  in 
(16)  after  N-1  terms.  We  find  the  N-1  coefficients  an  and  the  parameter 
y  by  collocation.  Thus  we  introduce  the  N  mesh  points 


)  **  0  (17) 


2?  dx  1  -x  ~  1  2  ^  d  ,  x 

TTcotgoe  —  +  —  e  cos0  cotga  cotga 


da'' 


ai =  1  =  (18) 

Using  (16)  and  (18)  we  obtain  [T(a)J  ,  [6(a)]  ,  [— ]  ,  [— ] 

o=o _  0=0  da  o=o  do  o=o 

2~  I  I  I  I 

d^Q 

and  [ — r)  terms  °f  the  coefficients  an  and  the  parameter  y . 

da  I 


Substituting  these  expressions  into  (16)  vre  obtain  N  nonlinear  equations  for 

the  N  unknowns  a  ,  n  ■  and  y.  We  solve  this  system  by 

n 

Newton's  method.  Once  this  system  is  solved  for  a  given  value  of  F,  the 
shape  of  the  jet  is  obtained  by  numerically  integrating  (6). 


i 

IV.  Discussion  of  the  result* 

We  use  the  numerical  scheme  described  in  Sec.  Ill  to  compute  solutions  of 

the  modified  problem  for  various  values  of  F  and  a.  The  coefficients  an  § 

were  found  to  decrease  rapidly  a  n  increases.  For  example 

a1Q  ~  10-2,  a2Q  ~  2.10-3,  a3Q  ~  3.10-*  for  F2  ■  0.02  and  a  “  10.  Most  of 

the  calculations  were  done  with  N  ■  50. 

In  Figure  4  we  present  values  of  y  versus  F  for  a  «  10.  As  F 
tends  to  infinity,  y  tends  to  zero.  As  F  approaches  zero,  y  oscillates 
infinitely  often  around  Figure  4  shows  that  there  exists  a  countably  j 

Infinite  number  of  values  of  F  for  which  y  -  The  solutions  corre¬ 
sponding  to  these  values  of  F  are  the  solutions  of  the  original  problem. 

Similar  results  were  found  for  other  values  of  a.  jj 

Vanden-Broeck  solved  the  problem  with  a  ■  (i.e.,  T  -  0).  He  found 

Y  »  V'  F  <  F  0.36 
2  c 

y  -  0,  F  >  F  ~  0.36. 

C 


(19) 

(20) 


Therefore  all  the  solutions  corresponding  to  F  <  Fc  a*e  solutions  of  the 
modified  problem  when  a  ■  *.  The  solution  defined  by  (19)  and  (20)  is  shown 
in  Fig.  4. 

The  flow  configuration  of  Fig.  1  can  also  serve  to  model  a  jet  emerging 
from  a  nozzle  whose  walls  consist  of  the  streamlines  IS  and  IS'  (see 
Vanden-Broeck  [2]).  Relation  (20)  shows  that  the  flow  leaves  the  wall  tan¬ 
gentially  for  o  »  •  and  F  >  Fc*  On  the  other  hand  Fig.  4  indicates  that 
the  flow  does  not  leave  the  walls  tangentially  for  a  *  10  unless  F  ■ 
Similar  discontinuities  in  slope  were  encountered  before  by  Vanden-Broeck 
[6,7]  in  his  analysis  of  the  effect  of  surface  tension  on  free-streamline 
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t 


Values  of  y  versus  F  for  a  ■  10.  The  broken  line  corresponds 
to  the  solution  defined  by  (19)  and  (20). 


flows.  The  recent  work  of  Vanden-Broeck  [10]  shows  that  these  unphysical 
discontinuities  could  be  removed  by  taking  into  account  the  thickness  of  the 
walls  and  the  finite  curvature  at  the  ends  of  the  walls. 

As  a  increases,  the  amplitudes  and  the  wavelengths  of  the  oscillations 
in  Figure  4  decrease.  In  Fig.  5  we  present  the  two  largest  values  of  F, 
which  correspond  to  solutions  of  the  original  problem,  as  functions  of  a  \ 
The  two  curves  approach  f  ■  F*  ~  0.23  as  a  tends  to  infinity. 

Our  numerical  results  suggest  the  following  general  result.  If  n 
denotes  an  arbitrary  positive  integer,  then  the  n**1  largest  value  of  F, 
which  corresponds  to  a  solution  of  the  original  problem,  approaches 
F  *  F*  «■  0.23  as  o  ♦  «.  Therefore  the  degenerancy  associated  with  T  ■  0 
is  removed  by  solving  the  problem  with  T  *  0  and  then  by  taking  the  limit  as 
T  ♦  0. 

These  results  imply  that  the  physically  relevant  solution,  when  T  *  0, 
is  not  the  solution  corresponding  to  F  -  Fff  ~  0.36  but  the  solution 
corresponding  to  F  -  F*  ~  0.23.  The  profile  corresponding  to  F  -  F*  and 
T  «•  0  is  shown  in  Figure  6. 

Collins3  found  the  experimental  value  F  ■  0.25.  In  addition  he  measured 
the  ratio  of  the  radius  of  curvature  at  the  top  of  the  bubble  to  the  width 
h  of  the  tube  and  obtained  the  value  0.305.  The  corresponding  ratio  for  the 
theoretical  profile  of  Fig.  6  is  0.32. 
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in  a  tube  of  width  h  is  considered.  The  acceleration  of  gravity  g  and 
the  surface  tension  T  are  included  in  the  dynamic  boundary  condition.  For 
T  ■  0,  the  results  of  Garabedian  and  Vanden-Broeck  show  that  a  solution 
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20.  ABSTRACT  (cont. ) 


exists  for  all  values  of  the  Froude  number  F  ■  U/(gh)^A  smaller  than  a 
critical  value  Fc  ~  0.36.  In  this  paper  accurate  numerical  solutions  with 
T  *  0  are  confuted  by  series  truncations.  The  results  show  that  for  each 
value  of  T  *  0,  there  exists  a  countably  infinite  number  of  solutions.  Each 
of  these  solutions  corresponds  to  a  different  value  of  F.  As  T  tends  to 
zero,  all  these  solutions  approach  a  unique  limiting  solution  characterized  by 
F  m  p*  ~  o.23.  Therefore,  the  degenerancy  associated  with  T  ■  0  is  removed 
by  including  the  effect  of  surface  tension.  In  addition  the  profile  corre¬ 
sponding  to  F  *  F*  is  found  to  be  in  good  agreement  with  experimental  data. 


